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Abstract We study a velocity-space quasilinear diffusion of the solar wind 
fS |' protons driven by oblique Alfven turbulence at proton kinetic scales. Turbulent 

fluctuation at these scales possess properties of kinetic Alfven waves (KAWs) 
that are efficient in Cherenkov resonant interactions. The proton diffusion pro- 
ceeds via Cherenkov kicks and forms a quasilinear plateau - nonthermal proton 
tail in the velocity distribution function (VDF). The tails extend in velocity space 
along the mean magnetic field from 1 to (1.5-3) Va, depending on the spectral 
break position, turbulence amplitude at the spectral break, and spectral slope 



D , after the break. The most favorable conditions for the tail generation occur in 



the regions where the proton thermal and Alfven velocities are about the same, 
Vtp/Va ~ 1- The estimated formation times are within 1-2 h for typical tails 
c/2 ' at 1 AU, which is much shorter than the solar wind expansion time. Our results 

r/2 , suggest that the nonthermal proton tails, observed in-situ at all heliocentric 

distances > 0.3 AU, are formed in the solar wind locally by the KAW turbulence. 



C/3 , We also suggest that the bump-on-tail features - proton beams, often seen in the 



proton VDFs, can be formed at a later evolution stage of the nonthermal tails 
by the time-of- flight effects. 
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Low collisionality of the solar wind (SW) and persistent activity of waves and 
turbulence make kinetic wave-particle interactions unavoidable in the SW mod- 
cling. The need for such kinetics has been emphasized in view of numerous 
non-thermal features observed in velocity distribution functions (VDFs) of the 
solar wind particles (sec Marsch (2006) and references therein). In particular, 
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nonthermal tails, beams, and temperature anisotropics are routinely observed in 
situ by satellites. 

Recent in-situ measurements have revealed that the SW turbulence at proton 
kinetic scales is dominated by Alfvenic waves (AWs) rather than the fast mode 
or whistler waves (|He et al.[ 120111 120121 |Podesta and Gary| I2UTTJ) . Furthermore, 
among Alfvcn waves, the most power, about 80%, was found in oblique quasi- 
perpendicular AWs, and the rest 20% in the quasi-parallel ion-cyclotron Alfvcn 
waves (ICAW) QHe et al., 2011 1. 



Properties of quasi-parallel ICAWs are understood much better than that of 
the oblique AWs at proton kinetic scales. ICAWs experience strong ion-cyclotron 
resonances and are considered as a main source for the perpendicular heating of 
the protons and heavier ions in the solar corona and solar wind (see reviews by 
|Hollweg and I senbcrg (2002]). [Marschl (|2006p and references therein). 

Properties of the dominant quasi-perpendicular AW fraction and its influence 
on particles VDFs are still unclear. Depending on the wavevector anisotropy 
k±/k\\ (perpendicular/parallel with respect to the mean magnetic field -Bo), less 
oblique AWs at proton kinetic scales are similar to ICAWs, while more oblique 
AWs are kinetic Alfven waves (KAWs). At intermediate k±/ku ~ Va/Vt p , AWs 
possess mixed properties of both ICAWs and KAWs and can be named ion- 
cyclotron KAWs - ICKAWs dVoitenko and Goossens, "20031 ). Since ICKAWs ex- 



perience both the ion-cyclotron resonance (typical for ICAWs) and the Chercnkov| 
resonance (typical for KAWs), they have more dissipation and generation chan- 
nels and can play an important role in the energy exchange between parallel and 
perpendicular degrees of freedom in the solar wind. 

The wavevector anisotropy at proton kinetic scales has been measured by 
Sahraoui et al. (2010). For perpendicular wave numbers up to k±p p ~ 2, the 
wave/proton-cyclotron frequency ratio appeared to be still low, w/fi p ~ 0.1, 
whereas the anisotropy is high, about k±/k\\ ~ 10, which is typical for classic 
KAWs. The dominant role of KAWs in the oblique AW fraction has been sup- 
ported further by two independent tests performed bv lSalem et al.\ (|2012p . who 
measured simultaneously two polarization ratios: compressional magnetic/total 
magnetic and total electric/total magnetic field perturbations. 

The interest in KAWs and related issues has risen recently in the context of 
turbulence dissipation in the solar wind and consequent heating of plasma species 
(Schckochihi n et ai\l\'2009^ ; Howes ( 201ip ; and references therein). IGoldreich and Sridhar[ 
(|1995p suggested that the MHD AW turbulence cascades along a critical balance 
path, which result in preferential generation of high perpendicular wavenum- 
bers, and hence KAWs ( |Schekochihin et al., 2009 ). An indirect confirmation of 



the critical balance was received recently by Gogoberidze, Chapman, and Hnat| 



(2012J), who proved that the critical balance leads to generation of the residual 
energy observed in the solar wind turbulence. 

A solar wind model elaborating Landau damping of KAWs has predicted the 
electron and proton heating rates ( |Howes, 201 1[ ), which agree well with empirical 
estimations bv lCranmer et al. (200SJ) at > 0.8 AU. At smaller radial distances 
the proton heating by Landau damping appeared to be insufficient and Howes 
(2011) suggested the proton cyclotron heating operating there. Because of the 
shorter incrtial range, the MHD turbulence at shorter radial distances ends up 
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with less anisotropic fluctuations at proton kinetic scales, which can produce the 
required proton heating via cyclotron resonant interaction. 

A widely used assumption of the M axwellian Landau damping can be violated 
in the solar wind by local deformations of VDFs of particles in resonant velocity 
ranges. |Borovsky and Gary| (|201ip have demonstrated that the electron Landau 
damping is still strong (although not Maxwellian any more) because the rela- 
tively high collisionality of the electron is capable of keeping their VDFs not far 
from Maxwellian. On the contrary, the less collisional protons experience much 
stronger departures from Maxwellian VDFs, which can significantly modify the 
proton Landau damping, or even cancel it. 

IChandran et qZ.1 (|2010p proposed another promising heating mechanism for 
the protons, stochastic acceleration by KAWs across Bq. This mechanism is not 
related to the kinetic wave-particle resonances, and is therefore not so sensitive 
to the local deformations of the proton VDFs. 

If the turbulence at proton kinetic scales consists of the ion-cyclotron and 
KAW fractions, as is suggested by recent observations, both these fractions 
should produce their own observed signatures - specific features in the proton 
VDFs. The most typical observed features of the proton VDFs are ( |Marsch, 1991[ ):[ 
(i) nonthermal tails or even secondary peaks along the magnetic field direction 
in all solar wind types (fast, slow, and intermediate), and (ii) total temperature 
anisotropics T p ± > T p » in the fast solar wind (with anisotropic proton cores), 
and Tpii > T p j_ in the slow (with isotropic cores) and intermediate (with both 
isotropic and anisotropic cores) solar winds. The parallel proton tails and beams 
are important contributors to the energy balance between parallel and perpendic- 
ular degrees of freedom. In the fast solar wind they reduce the total temperature 
anisotropy to T p ±/T p \\ > 1 (and at times even reverse to T p || > T p ±) despite the 
strongly anisotropic cores of the proton VDFs, T££ e /T£? ie = 2 — 3 (see Figure 
8.4 by Marsch, 1991). In the intermediate and slow solar wind they dominate the 
energy balance making T p \\/T p ± > 1. Among them, the nonthermal tails prevail 
QMarsch et al., 1982[ ). The secondary peaks (bumps-on-tails) are less frequent 
and develop on the tail background mainly in the fast solar wind. Therefore 
the bump-on-tail seems to be a later evolution stage of the proton VDF as 
compared to the tail. At the same time, persistent nonthermal tails without 
bumps are observed at all distances, which suggest that the bumps are not 
necessary developing on the tails within the solar wind expansion time-scale; 
their development is therefore somehow constrained or require more time. 

The anisotropic cores of the proton VDFs with TfJ™ > T c n re have been 
associated with the ion-cyclotron heating by ion-cyclotron resonant waves, and 
ICAWs interaction with solar wind protons and heavier ions have been stud- 
ied extensively QHollweg and Isenb erg (2002); Marsch (2006); and references 
therein). It was in particular shown that the ion-cyclotron pitch-angle diffu- 
sion can be responsible for the shaping and perpendicular heating of the core 



proton VDFs (Marsch and Tu, 2001), for the cross-field heating of heavier ions 



(|Galinsky and Shevchenko, 2012 1 , and also for the pitch-angle scattering of the 



nonthermal proton tails (Marsch and Bourouaine, 20111 



Wave mechanisms producing nonthermal proton tails and bumps-on-tails 
(beams) in the solar wind received much less attention. The first one, proposed 
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by |Tu, Wang, and Marsch| (|2002|) . relied upon the proton cyclotron-resonant 
diffusion driven by a specific ion-cyclotron mode supported by helium ions. 
Another mechanism, nonlinear proton trapping in the parallel KAW poten- 
tials, has been suggested by I Voitenko and Goossensl(|2006p . Later on, analytical 
( jPierrard and Voitenko, 2010[ ) and numerical ( |Li et al., 2010| |Osmane, Hamza, and Meziane, 2010[ )[ 
studies have confirmed that a KAW with reasonable amplitude can produce 
observed proton beams in the solar wind conditions. 

Proton trapping and beam generation can also be produced by electrostatic 
waves carrying parallel electric fields. Recent simulations of the nonlinear decay 
of circularly polarized Alfven waves have demonstrated that the product electro- 
static waves can trap and accelerate protons along Bq, thereby creating a beam 
(Araneda, Marsch, and Vinas, 2008; Matteini et al., 2011||Valentini, Perrone, and Veltri, 2011[).[ 



Ion-acoustic waves ( Araneda, Marsc h, and Vinas, 2 008 ; Mattei ni et al., 2011 1 and| 



"ion-bulk" waves (Valentini, Perrone, and Veltri, 2011) have been discussed as 



mediators transferring energy from AWs to plasma particles. Alternatively. IRudakov et qZ. [ 
(|2012p used a quasilincar approach to study the proton diffusion driven by 
a strongly dispersive KAW spectrum. This process can produce a quasilinear 
plateau in the form of nonthermal tail in the proton VDF. 

Earlier kinetic models of the SW electrons, accounting for the influence of ex- 



ternal macroscopic forces and Coulomb collisions ( Pierrard, Maksimovic, and Lemaire, 1999] [ 



Lie-Svendsen and Leer, 2000 Vocks and Mann, 2009), have been recently im 



proved by the inclusion of terms due to whistler turbulence (Pierrard, Lazar, and Schlickciser, 2 011[ ).[ 

It was concluded that the whistler turbulence, if exists in the solar wind, could 

lead to the velocity-space diffusion of electrons and form nonthermal tails in the 

electron VDF. Another source for the velocity-space electron diffusion has been 

studied bv IRudakov et all (|2011[) . who assumed a spectrum of kinetic Alfven 

waves (KAWs) instead of whistlers. Moreover . IRudakov et aLI(|201ip argued that 

the nonlinear KAW scattering off electrons is more efficient than the nonlinear 

wave- wave interaction in the short- wavelength KAW range. 

In the present paper we study the influence of the observed kinetic-scale 
turbulence on the evolution of the proton VDF. First we derive the Fokker- 
Planck diffusion terms for protons under the action of Alfvcnic turbulence. We 
are interested here in the influence of the KAW component, which is dominant 
in the kinetic-scale SW turbulence, on the evolution of the proton VDF. The 
Fokker-Planck diffusion coefficients are presented in terms of quantities measured 
in-situ in the solar wind. In particular we refer to magnetic field measure- 
ments that provided spectral indexes of the SW turbulence at proton kinetic 
scales ( |Alexandrova et al., 2009| |Sahraoui et al., 2010[ ), MHD/kinetic spectral 
break wavenumber, and the turbulence amplitude at the break wavenumbcr 
( |Smith et al., 20 06; Markovskii, Vasquez, and Smith, 2008). Analytical and nu- 
merical estimations are given for the evolution of nonthermal tails in the proton 
VDFs. 
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2. Wave Model 

The classic KAW dispersion can be presented as 

u k = k {l V A K, (1) 

where k is the wave vector, || (J-)) mean components parallel(perpendicular)to 
Bo, K = K (fi) is the KAW dispersion function, and /j = p p k± is the dimension- 
less perpendicular wavenumber. 

The kinetic dispersion function for KAWs was derived by|Hascgawa and Chen| 
dM5D: 



where A = A (fi) = I (/i 2 ) exp (— /i 2 ), I (fJ. 2 ) is the zero-order modified 
Bcsscl function, and T G ( p ) is the electron(proton) temperature. This expression 
was obtained in the following wave and plasma parameter ranges: k 2 jk\ <C 1, 



,2 



2 



u>l/nl « i, v 2 p « vl < (« fc /fcu) « v Tc 

The KAW dispersion function gives the KAW phase velocity, Vk = ^>k/k\\, in 
units of Alfven velocity: Vk/V A = K (fi). In the MHD limit /i-)0 KAWs become 
dispersiveless Alfven waves, Vk — > Va, but with growing /z the KAW phase 
velocity deviates from the Alfven velocity significantly. The inverse function 
K^ 1 (14/Va), which we will need later in the diffusion coefficient, is impossible 
to find from Equation ([2]) analytically in general case. Explicit expressions for 
K^ 1 (V^/Va) can only be found in two asymptotic regimes of weak ([i 2 <C 1) 
and strong (/i 2 3> 1) wave dispersion. 

To simplify the problem, one can use a Padc approximation for Aq (/«), 



f 

1+"/! 2 



A o(m)-T-^^> ( 3 ) 



resulting in the following KAW dispersion function 
and its derivative 




T n 



li 



dKp/d ^{ 1 + ^l)W (5) 

These expressions can also be found in the framework of two-fluid MHD 
plasma model and provide a good approximation for KAWs in the whole range 
of fi if the plasma beta /? = (l + T ||/T p j_) V^ p /Vl <C 1. As is seen from Figure 
1, the difference between the kinetic KAW dispersion K (fi) (solid line) and its' 
Pade approximation Kp (fi) (dashed line) is insignificant for T e u/T p j_ at 1. 

Finite-/? effects come into play in the solar wind at 1 AU, where the typical 
value |8 « 1 (this corresponds to the plasma/magnetic pressure ratio «1). The 
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K(n) 




Figure 1. Comparison of three models for the KAW dispersion function K (p.): (1) kinetic 
dispersion by Hasegawa and Chen (1976) for /3<1 (solid line); (2) Pade approximation, that 
corresponds to the two-fluid MHD dispersion for j3 <C 1 (dash line); (3) two-fluid dispersion 
with finite-/? effects for /3ii = /3j_ = 0.5 (dot line) and its asymptote (bottom solid line). In all 
cases T e ||/T p x = 1. 



KAW dispersion function Kp, accounting for the finite- j3 effects of magnetic and 
plasma compressibility flVoitenko and Goossens, 2002[ ), can be modified for the 
case of anisotropic temperatures: 



K} 



1 /3|| , 1+Teii/TVL 2 

1 + l+/3x + l+0± V 



\ 



l+/3_i 



1 + T e\]/ T pX 2 
1+0X ' 



^11 



l + /3j 



(6) 



A . This dispersion is also shown in 



where /3||, ± = {l + T clu /T pl± )V^JVl 

Figure f (dot line). The corresponding KAW wavenumber [i can be expressed 

via Vk explicitly, 



M = 



l + /3j_ (&) - 


-d+ fh \ ( Vk Y i /s " 


. 1 + T e ,|/T p± 


(t) ! 



(7) 



At large /x, the KAW dispersion (phase velocity) can be decreased significantly 
by the finite-/? effects. On the contrary, in the weakly dispersive range fi < 1 
such moderate values as j3 = 0.5 do not change the KAW dispersion significantly 
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and can be approximated there as 



, i (Vk 

Va 



A* = k p I — 



fci'-' 



\ i 






(8) 



3. Proton Diffusion in the KAW Turbulence 

The parallel component of the KAW electric field E z \^ (k = (k x ,k y ,k z ) is the 
wave vector, z || Bo), makes KAWs efficient in Cherenkov interaction with 
plasma particles. Following ( |Voitenko, 1998| [Voitenko and Goossens, 2003| ), we 
account for the proton VDFs modifications under the action of KAWs in the 
framework of quasi-linear theory. We consider an axially symmetric (with respect 
to the background magnetic field B ) problem, where both the KAW spectrum 
and the particle velocity distributions are independent of respective polar angles 
in the cross-field plane. 

In general, the proton distribution function F = F(r, V, t) in the solar wind 
obeys the Vlasov collisional equation 

9 ,r 9 9\ n fdF\ fdF\ 

where r and V are respectively the position and velocity vectors of the protons, 
t is the time, and a is the proton acceleration under the action of the exter- 
nal forces: the macroscopic electric force ZeEo, the gravitational force mg and 
the Lorentz force oc v x Bo- The right-hand side represents the velocity-space 
proton diffusion due to Coulomb collisions, (dF/dt) c , and due to wave-particle 
collisions, 

dF\ ^^_ DK dF_ 
dt J A dV z dV z 

In the diffusion coefficient D A we integrate over the polar angles in the cross- 
Bo plain, implying that |-E z k| is axially symmetric, but we keep the dependence 
on the perpendicular component of the particle velocity: 

2 

D A = ^Y, 6 ( UJk - k * V ^ J o\ E ^\ 2 ( 10 ) 

P k 

where the zero-order Bessel function J = J (k±PpV±/Vr p ) reflects the fact that 
the actual electric field of KAWs acting on the protons is reduced because of 
the averaging over their cyclotron orbits with the velocity-dependent gyroradius 
Vj_/f2 p = p p V±/Vt p - Other notations are as follows: k z and k± are the wave 
vector components parallel and perpendicular to Bq, k\ = k^. + ky, p p = Vt p /Q, p 

is the proton gyroradius, Vr P = \jT p jm p thermal velocity, T p temperature, 
m p proton mass, q p proton charge, and J^ k — 2tt J dk z J dk±k±. Dirac's delta- 
function S (uik — k z V z ) follows from the resonant character of the Cherenkov 
wave-particle interaction. 
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The simplest solar wind models consider only the effects of the external forces. 
In such models called exospheric, the right-hand side interaction terms in Equa- 
tion ([9]) are neglected so that an analytic solution of the equation can be obtained 
( |Lamy et al.\ (|2003p and reference therein) . Coulomb collisions make equations 
too complex and require numerical simulations. Expressions for the Coulomb 
diffusion term in the Vlasov equation for protons and results of numerical simula- 
tions accounting for Coulomb collisions are presented in ( |Pierrard and Voitenko, 2012[ ).[ 

Starting from Equation (|10[) , the wave-particle diffusion term accounting for 
turbulence properties is derived in next sections. We especially pay attention to a 
non-monotonous dependence of the diffusion coefficient on the parallel velocity, 
which follows from the non-monotonous dependence of the spectra of parallel 
electric fields on the KAW perpendicular wavenumber. 

3.1. Fokkcr-Planck Diffusion Coefficient in Terms of KAWs' Magnetic Fields 

Since the parallel electric fields E z \^ in KAWs are relatively weak and difficult 
to measure, it is instructive to express D A in terms of KAWs' magnetic fields. 
To this end one can use the KAW polarization relation 

B ± T p _l c k±K' [ ' 

Then the 2D velocity-space diffusion coefficient for the protons can be reduced 
to the following integral in the normalized wavenumber space: 

T ™ n t/2 / ... e /. , „-i\ ^4 !dw\B tlu \ 2 



D A = ^^-SlpViJdfiSifi-K 



p± p b J ^ v > K^dK/dfx B 2 



Jp_L 



^4 jdvv\B, 



I 2 



K 2 dK/dn Bl 



2 



(12) 



=K~ 1 



The dimensionless parallel wavenumber v = d p k z (6 P is the ion inertial length) is 
introduced in Equation (|T2"j) . and the following relation for the Dirac S- function 
is used: 

5{K -V z /Va) : 



where K~ l = K~ x (V z /Va) is the inverse A'-function of V z /Va and Vs 



2 



^/T ez /m p is the ion-sound speed. The normalized spectral density, \B^ V 

> jllJ 



|Bj_k| , is defined such that J dk z J dk±k± \B±\ t \ = / dv J d^ii\B^ 



Note that \B^\ has the same dimension as Bq. 

Now we have to express the integral L dvv \B^ V \ in terms of the integral 
power at //, \B^\ = J dv \B^ V \ . This last quantity is related to the unidirectional 
spectral density W M measured by spacecraft: /i \B^\ = W^. The perpendicular 
wavenumber fi is related to the measured spacecraft-frame frequency / through 
fism9vB — (Vtp/Vsw) (///p), &vb is the angle between the solar wind velocity 
V S w and B . 
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In principle, an unknown spectrum of parallel wavenumbers v contributes to 
the turbulence level at every particular [i. To simplify the problem, we take 
into account the following two facts. First, theory and observations suggest that 
the turbulence cascade proceeds along a path in wavenumber space which is 
defined by the "critical balance" between the linear and nonlinear time scales 



( Goldreich and Sridhar, 1995 1. In accordance to the critical balance condition, 
every particular /x has its own v^ where most of the spectral density is con- 
centrated: \B^ V \ = \B^\ 8{v— v^,), where |B M | is the turbulence level at the 
perpendicular wavenumber (i. Second, the turbulence fluctuations are anisotropic 
with kj_ 3> k z . Denoting the anisotropy factor k z /kj_ = a (/cj_), the critical bal- 
ance reads v^ = y^a (/i) \i. Hence we estimate J °° dvv \B^ V \ ~ y~a (/i) /i \B^\ 



^a(/i)iy M ,and 



n A 2 



T f 



\ 3/2 



where /iv = A" 1 (V Z /V A ). 



V S V A 



:(//)/X 



3 t2 
J 



w„. 



K 2 dK/dn B'l 



(13) 



M=/iv 



3.2. Diffusion Coefficient for Double-kink Power-law Turbulence Spectra 



Further progress is possible if we know a particular shape of the turbulent 
spectrum of the Alfvenic turbulence. Numerous in-situ observations indicate 
that the power law spectra Wy, oc \i~ v are typical in the solar wind, but the 
spectral index p is different in different ranges of fi. In the MHD turbulence 
range /i < \x\ the power law index p = p\ ~ 1.7, close to the Kolmogorov 
value. In the intermediate wcakly/mildy dispersive KAW range \x\ < p, < 1 
the spectra are much steeper, with p = pi varying from p 2 — 2 to p 2 = 4 
(Leamon et al. 1999; Smith et al. 2011; Sahraoui et al, 2010). In the strongly 
dispersive KAW range, fi ^> 1, the power index approaches the value p = ps ~ 
2.8 QAlexandrova et al, 2009| |Sahraoui et al, 2010] ). All the above double- kink 
behavior of p (/x v ) can be modeled by the following piecewise function: 

P = p{fJ>) = Pi + (P2 - Pi) H (fi - m) + (p 3 - p 2 ) H (fi - fi 2 ) , (14) 

where H (x) is the Heaviside step function, H (x) = for x < and H (x) = 1 for 
x > 0. The spectral break (more precisely, the first spectral kink) at [i\ = 0.1—0.5 
separates MHD and (weakly/ dispersive) kinetic ranges of Alfvenic turbulence. 
The second spectral kink, that occurs at \x 2 = 2 — 3 (see spectra measured by 
ISahraoui et al.\ (|2010|) . separates weakly and strongly dispersive KAW ranges 



( Voitenko and De Keyser, 2011 ). 

With Equation (|14[) . the measured spectral density can be written in the 
following general form: 



Wy = Wy 



>® 



(P3-P2)H(fl-fl 2 ) 



-(Pl+(P2-Pl)H(n-y 1 ) + (p3-p 2 )H(fJ,-IJ,2)) 



(15) 
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where W^i is the turbulence spectral density at /i = /Ji (note that W^\ has 
dimension of Bq). The re- normalization of the spectrum (Equation (|15[) ^ at \i > 
/X2 ensures it is continuous through the second kink at /J, = /«2- 

With Equation p3]) the velocity diffusion coefficient (Equation ([T5|l ) attains 
the following general form: 



D A = ir 2 n p V s V A 



(fc) 



3/2 



1 + ^ 



,2^„ UMl '«-»)^-«)/ fc \ 2 *) 



Va 



v^ " -"^-. ( l6 ) 



n 



where /z v = p-(Vk = Vfi) is the resonant KAW wavenumber defined by Equation 
and p (jUv) is given by Equation (p~4|) . In what follows we will put a (p.) ~ O.l 
at ion kinetic scales /j«l ( Sahraoui e£ a/., 20I0[ ), and neglect its relatively slow 
variation in the resonant wavenumber range. 

The diffusion due to KAWs is strongest for the protons with velocities around 
V z c± Vi > Va, which are Cherenkov- resonant with the spectral peak of the KAW 
parallel electric fields. In the plane _L Bq, the diffusion is maximized for the core 
protons at V± < Vt p , where the reduction due to proton-cyclotron gyration is 
minimized. In the -Bn-parailcl direction, the diffusion coefficient decreases very 
fast as V z —)■ Va, and less fast with V z growing in a more extended velocity 
range from V\ to several Alfven velocities. These features are compatible with 
the physical picture of the Cherenkov resonant interaction between gyrating 
particles and waves with finite cross- Bq length scales. 

Since there are no positive /i = K^ 1 (V z /Va) for V z /Va < I, D A = there 
(this fact reflects the absence of resonant KAWs for particles moving with sub- 
Alfvcnic velocities V z /Va < 1). One should note, however, that even sub-Alfvcnic 
protons can be affected by the turbulence fluctuations via two effects: (i) nonlin- 
ear broadening of the Cherenkov resonances, and (ii) increased Coulomb diffusion 
of the protons at the steep VDF slope at V z < Va- The minimum velocity of the 
affected protons resulting from the former effect is analytically estimated below 
in Equation (|22)l . 

It is interesting to note that the spectral slope p (/x v ) < 2 in (fbT| can make 
D a still increasing function of /i v in the dissipation range between {i\ and 
fi2- This follows from the increase of the parallel electric fields in KAWs that 
is faster than the power-law decrease of the turbulence amplitudes. Therefore, 
depending on p(fi v ) ;=s 2, the diffusion coefficient attains a maximum value at 
/i v = pi (with p(/i v ) > 2) or at /i v < /12 (with p{p. v ) < 2). Such behavior of 
the diffusion coefficient could have an important consequence that the less steep 
dissipation-range spectra result in longer nonthermal tails in the proton VDFs. 
However this could only be true for the cases where the turbulence amplitude is 
the same at p, = p,\. 

In fact, as was found by ISmith et al.\ (|2006p in the inertial range, the larger 
spectral fluxes e are followed by the steeper dissipation-range spectra P2, such 
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Figure 2. Reduced diffusion coefficient as function of normalized parallel velocity V = V z /V\ 
for p2 = 2 (dot line), P2 = 3 (dash line), and p2 = 4 (solid line). Turbulence level at the first 
spectral kink is adjusted for each p2 using the scaling Wul °c p 7 ' 4 . 



that p2 = 1.05e 009 . Applying this relation at the end of inertial range, we find 
the scaling 

W^txpl 4 , (17) 

relating the spectral power W^\ at the break \i = [i\ and the spectral index pi 
above it. If the spectral powers in two cases (1) and (2) arc W ^ and W^ x , 

then the ratio WAi /W^i = (^2 lv\ J ■ Elaborating this scaling relation, a 
reduced ID diffusion coefficient is shown in Figure 2 for three spectral indices 
P2 = 4, 3, and 2 from top to bottom. The reduction was performed by the 
averaging of Jq in (|Tt5)) in the cross-field velocity plane. With the correlation 
found by Smith et al. (2006), the diffusion coefficient, in average, is larger for 
steeper dissipation-range spectra, even if they decrease with V z faster. 

The described above properties of the diffusion coefficient indicate that ki- 
netic Alfven turbulence can be efficient in diffusing the protons in the velocity 
range covered by nonthermal tails of the proton VDFs. This is supported by the 
analytical and numerical estimations presented below. 



4. Evolution of Nonthermal Tails in the Proton VDFs 

We investigate here what nonthermal proton tails can be formed by their quasi- 
linear diffusion along the kinetic turbulent spectrum described above. Let us 
consider a local evolution problem in the solar wind neglecting external forces, 
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spatial variations, and Coulomb collisions in Equation ©. To estimate the for- 
mation time for a particular tail, we have to know the tail extension in the 
field- aligned direction Vmax- Observations show that the tail extension Vmax 
rarely exceeds 3Va, which allows us to simplify the problem accounting only 
for the intermediate KAW spectrum in the diffusion coefficient (Equation (TT51) ) . 
The cross-field velocity spread in the tails is approximately the same as in the 
proton core, ~ Vt p j_ (within the multiplier of the order 1). 

Accounting for above, we reduce Equation ([9]) to the ID diffusion equation 

m ~ W D W Z (18) 

for the reduced ID distribution function / = f (V z ) = 2n J dV±V±F (V Z ,V±). 
The ID diffusion coefficient in the range Va < V z < 3Va, 

p A (l + T e2 /T p± ) ^ V z Bl' K ' 

follows from Equation fj 16f) after averaging over Maxwcllian distribution in the 
cross-field plane: I Jq (mv^) ) = A (fJ-l)- 

For small /ii the corresponding resonant velocity V\ = VaK (fii) is close to 
Va and the remaining velocity interval Va < V z < "V\ is narrow. The detailed 
behavior of the diffusion coefficient D is unimportant there because D decreases 
fast to with V z — > Va, and the nonlinear resonance broadening of KAWs with 
/i ~ hi is more important. 

4.1. Quasilinear Plateau Model 

In this section we study the quasilinear plateau evolution under the influence of 
the KAW turbulence. To this end we start from the initially isotropic Maxwcllian 
VDF for the protons 

1 f V? \ 



ft=o = hi = r— exp . . . 

v^rl/ T p|| \ ZK x P || / 

The diffusion coefficient (Equation (fTO)) ) and hence the velocity-space diffusion 
are very non- uniform, as is seen from Figure 2. The diffusion is fastest in the 
vicinity of Va , where the diffusion coefficient is maximal, but decreases about two 
orders at 3Va- The quasilinear saturation in this situation occurs first around 
this maximum, where an initial local plateau is formed as is shown in Figure 3 
(dashed line). The resulting ID proton VDF can be modeled by the piece- wise 
function with the local plateau 

/ = /pi = v fM v (20) 

'max ^min 

in the velocity range Vmin < V z < V max , and / remaining Maxwellian, / = /m, 
outside the plateau. 
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Figure 3. Model proton VDF with a step-like quasilinear plateau extending to V ma x = 2.5Va 
at the late evolutional stage (an early evolutional stage is shown by dash line). Numerical 
solution of the diffusion equation (Equation JTH)) for the same time is shown by stars. A good 
correspondence is observed between the model and numerically calculated VDFs. Details on 
the plasma and turbulence parameters used are given in the text. 

The quasilinear saturation spreads in time to higher velocities, shifting the 
plateau front V ma x = Vmax (£) forward to higher velocities as is also shown in 
Figure 3 (solid line). The normalized proton number density in the plateau can 
be expressed in terms of error function Erf: 



V 



"-pi 
n 



dV z f M = 



Erf Gfe 



-Erf 



</2V T 



(21) 



It increases with time, but the plateau height decreases. 

It is important to note that the minimum velocity of the protons involved 
in the plateau formation, VJ n i n , can be less than the minimum phase velocity of 
resonant KAWs, Va- As follows from the nonlinear broadening of the Cherenkov 
resonance for a KAW with wavenumber kj_ and amplitude Bk, the minimum 
velocity of the affected protons is 



V^ 



V A 



= ji + kip^- 




T «\\/ T P 



k±PT Bk 



V A ^Tp^/Tpii+Teii/TpiiY y/l + k\p* T B 



(22) 



This formula follows from the condition that the proton kinetic energy in the 
wave frame is equal to the KAW potential barrier. For T c \\/T p \\ = T p ±/T p \\ = 1, 
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and the typical KAW amplitude Bk/Bo = 0.05 in the vicinity of the first spectral 
break, at fcj_PT = yl + T c \\/T p ±k±p p ~ 0.3, the minimum velocity reduces to 
Vmhi ~ 0.8Va- The modified V m [ n < Va can significantly increase the proton 
number density in the tail, and hence the energy in the parallel degree of freedom 
in comparison to the case V m i n = Va • On the other hand, such decrease of V m { n 
does not influence much the evolution time scales of the tails. 

4.2. Evolution Time Scale of the Quasilinear Plateau 

The time evolution of the plateau spreading can be found from the diffusion 
equation (Equation (0) as follows. First, we note that the time derivative of 
/pi can be expressed as 

d/pi _ /pi — /m (Vmax) dV max . . 

dt ~ V max -V mia dt l ' 

Using expression (|23|) in the lhs of the diffusion equation (Equation (fTg|) ) and 
then integrating it over the plateau velocity range, from V m i n + to V" ma x + 0, we 
obtain a first-order ordinary differential equation for Vmax as function of time t: 

ir p /t/ \\ "'max r, rt/ \ & JM {Vmax.) 
(Jpl — /M i^maxJJ — -T — L) ( V max ) — . 

dt o v max 

This equation can be solved easily by separating variables and integrating from 
t = to t: 

. _ v-2 f m " X /pl (Xnax) //m (Knax) — ^ jt// CO/l\ 

^ - ^TpH / y, D( y, v d^max- (^j 

Let us estimate the plateau formation time at 1 AU, where £7 p ~ 1 s _1 , 
T cz /T p± ~ T cz /T pz ~ 1 and Va =^ V Tp ~ 50 km s" 1 . With the turbulence 
amplitude \B±\ /Bq ~ 0.05, spectral break wavenumber /hi ~ 0.6, and spectral 
slope p2 = 3 just above the break, (|2"4l gives the time about half an hour for 
formation of quasilinear plateau with the length equal to the average tail length 
V max = 1.75Va- Since this time is much (more than 2 orders) shorter than the 
solar wind expansion time igw (2-4 days at 1 AU), such tails can be easily 
generated by the turbulence locally in the solar wind. 

A long tail with Vn iax reaching 2.5Va is developed by the P2 — 2 turbulence 
at t = 14 h, when the local approximation is still marginally applicable. The 
corresponding model solution is shown in Figure 3 by the solid line. Further- 
more, to verify the model, we solve numerically the proton diffusion equation 
(Equation (fTS]) ) with the zero-gradient Neumann boundary conditions at the 
lower (V z /Va = 1) and upper (V z /Va = 5) boundaries. The value for the upper 
boundary is chosen in such a way that its variations does not affect solution in the 
velocity range V z /Va = 1 — 3. This solution is also shown in Figure 3 by stars. A 
good correspondence between numerical and analytical solutions indicates that 
the model (Equations (|2T)|) - (j2"4"|) ) is sufficient for quantitative estimations of the 
quasilinear plateau height, length, and evolution time-scale. 
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For longest tails with V max ~ 3Va, the tail formation time approaches the 
solar wind expansion time isw = 2 — 4 days. This means that the local analysis 
for such tails is inapplicable, but does not mean that such tails cannot be 
formed by the turbulence. To study such tails one needs to solve a nonlocal 
problem incorporating radial dependencies of the solar wind plasma and turbu- 
lence characteristics (this work is in progress). Again, even for such long tails as 
Vmax — 3Va, the problem can be made local with slightly increased turbulence 
amplitudes and/or spectral break wavenumbers, say to \B±\/Bq ~ 0.07 and 
lii ~ 0.8. Such variations are within the ranges of measured values (see Figure 
2 by Markovskii, Vas quez, and Smith| (|2008l . In addition, as we will see below, 
the decrease of the spectral slope P2 (with fixed amplitude B\) has the same 
effect. 

Since the plateau formation time given by Equation (|24|) is quite complex func- 
tion of plasma and turbulence parameters, its radial dependence require further 
investigations using particular plasma and turbulence models. For example, with 
approximately constant break frequency (constant fcj_i), the plateau formation 
time decreases with decreasing R as, roughly, ~ Va ~ R, the same as isw- In 
this case the locality condition can be still satisfied much closer to the Sun, 
where Va is larger and hence longer tails can be generated. Several numerical 
estimations with reasonable plasma and turbulence parameters have shown that 
the tails Vm ax ps 2Va can be easily generated, within half an hour, at such short 
heliocentric distances as R « 0.1 AU, where isw « 5 — 7 h. The conditions are 
even more favorable with constant fi\. These estimations are compatible with 
Helios observations of nonthermal proton tails at all distances from 0.3 to 1 AU. 

Figure 4 shows the time evolution of the plateau boundary V max =V rnax (t) in a 
turbulence with the spectral power W/i ~ 0.2 nT 2 Hz -1 at the break frequency 
/i ~ 0.3 Hz. These values are well within ranges of measured values and close to 
the average ones (see Figures 1 and 2 by |Markovskii, Vasquez, and Smith|(|2008p . 
Five curves in Figure 4 correspond to five different spectral slopes P2 =2, 2.5, 
3, 3.5, and 4 (from top to bottom). Longer dashes indicate larger p2, except 
for the solid curve with p2 =3. It is seen that the formation time t = £ql for 
the plateau V max = (1.5— 2.5) Va is considerably shorter than the solar wind 
expansion time, £ql/£sw = 0.01 — 0.3 (we assume an average isw = 2.7 x 
10 s). With these values, covering typical proton nonthermal tails observed in 
the solar wind, the quasilinear diffusion process can be considered in the local 
approximation and our results are self-consistent. As the tail length approaches 
3Va, the corresponding evolution time scale approaches isw for the shallowest 
spectra. This suggests that in average plasma conditions the tail length can 
hardly exceed 3Va, except for the cases of high turbulence levels combined with 
shallow spectra. 

In Figure 5 we are looking for the favorable conditions facilitating local 
production of nonthermal proton tails in the solar wind. To this end we fix 
the plateau front at V max = 1.5 and plot the corresponding formation time 
as function of Vt p /Va (which is proportional to the square root of the proton 
plasma beta). The spectral indexes are P2 =2.5, 3, 3.5, and 4 for the four curves 
from bottom to top. The most favorable Vt p /Va for each curve is that one 
at which i/isw attains minimum. The shallowest spectra appeared to be most 
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Figure 4. Normalized tail length Vmax/VA as function of normalized time i/isw- Spectral 
slopes are p2 = 2. 2.5, 3, 3.5, and 4 from the top curve to bottom. All other plasma and 
turbulence parameters are the same for all curves (see text). 



efficient at Vt p /Va — 1-5, whereas the steepest spectra at Vt v /Va < 1- These 
values are not much different for all spectral slopes of interest, so that we can 
take the optimal value Vt p /Va ~ 1, which is close to the average value in the 
solar wind at 1 AU. This means that the turbulence at 1 AU become more 
efficient in producing nonthermal tails. 

In general, as is seen from Figures 4 and 5, the shallower kinetic spectra are 
more efficient in the tail production provided all other turbulence characteristics 
are fixed. However, not only the spectral slope, but also the level of turbulence 
is highly variable in the solar wind. Moreover, as was noticed by ISmith et al.\ 
(|2006p , there is a positive correlation between the spectral flux before the spectral 
break and the spectral slope after it. The corresponding scaling in terms of spec- 
tral power can be written as Wp\ oc p\ . Rescaling the turbulence amplitudes 
corresponding to P2 =2, 2.5, 3.5, and 4 slopes, we plot the resulting plateau front 
velocities for each p 2 in Figure 6. Now we obtain an opposite ordering of V max 
with p 2 , indicating that longer tails are generated by steeper spectra. Another 
new feature is less scattered final values of V max as compared to the case of 
uncorrelated amplitudes and slopes. 

4.3. Numerical Example 

At present, there are no data on the proton VDF and kinetic-scale turbulence 
measured simultaneously. The problem with the Cluster spacecraft is that they 



NUMREF: proton . turb . diffusion; 5 April 2013; 2:12; p. 16 



Velocity-Space Proton Diffusion in the Solar Wind Turbulence 



1.00 



t 



0.50 - 



0.20 - 



0.10 



0.0 



0.0 




0.01 - 



10.0 



Figure 5. Formation time of the tail with V m ax = 1-5Va as function of the proton ther- 
mal/Alfven velocity ratio Vrp/^A f° r th e spectral slopes p2 = 2.5, 3, 3.5, and 4 (from top to 
bottom). Other parameters are as in Figures 3 and 4. 



are capable of measuring the turbulence at sufficient frequency resolution, but 
the proton velocity-space resolution is insufficient in the solar wind conditions. It 
is nevertheless interesting to estimate a nonthermal proton tail that could result 
from a turbulence measured by Cluster. 

In a particular case of Cluster measurements 10 January 2004 from 06:15 to 
06:25 UT, reported by Sahraoui et al. (2010), the background plasma parameters 
are: B = 10.2 nT, n = 16 cm" 3 , T c = 10.4 eV, T p = 31 eV, V sw =548 
km s , angle between the solar wind velocity and magnetic field 9bv = 67°. 
Therefore, the proton thermal velocity in this case is about the same as the 
Alfven velocity, Vt p — Va — 55 km s _1 , proton cyclotron frequency f2 p ~ 0.94 
s _1 , proton gyroradius p p ~ 59 km, and solar wind expansion time isw = 76 
hours. The key turbulence characteristics are: spectral index at the proton kinetic 
scales p = P2 — 4, spectral break \i\ = 0.4, and the turbulence spectral density 
W^i ~ 0.15 nT 2 at /x = Mi- 

The dependence V max (t) for this case is given by the bottom line in Figure 4. 
One can observe that the turbulence should generate a nonthermal proton tail 
with V max ;$ 2Va at 1 AU. The proton VDF could not be measured by Cluster, 
but the average parallel energy of the protons was found to be larger than the 
perpendicular one, producing the effective temperature anisotropy T p \\/T p ± < 
1.5. This anisotropy cannot be produced by the quasilinear evolution of the 
initially isotropic Maxwellian VDF in the velocity range from Va to 2Va- To 
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Figure 6. Normalized tail length V max /VA as function of normalized time t/tsw f° r the same 
parameters as in Figure 4 except for the turbulence amplitudes are now adjusted by the scaling 
M^ul °c v\ ' 4 - in this case the order of curves is reversed, P2 = 2, 2.5, 3, 3.5, and 4 from bottom 
to top, which means the longer tails are generated by steeper spectra. 



obtain such anisotropy, one needs either more protons in the tail, or longer tail 
(Vmax > 2Va), or both. The additional protons can be delivered to the tail by 
the nonlinear broadening of the Cherenkov resonance at /x ^ fj,\, as is given by 
Equation (|22p , and/or by the enhanced collisional diffusion at the steep velocity- 
space gradients at V z ~ V m - m , where VDF undergoes a sharp transition from the 
Maxwellial VDF to the plateau-like. Most probably, both these processes have 
non-negligible effects on the tail density, and act synergetically. Again, we would 
like to stress that the tail length is not much affected by these processes. The 
longer tail could only be generated at shorter radial distances to the Sun, in 
more favorable conditions. 



5. Conclusions and Discussion 

We studied a Fokker-Planck proton diffusion in the solar wind induced by 
Alfvenic turbulence observed recently at proton kinetic scales (jSahraoui et al.{ 
mm IHe et al\ IMll20l2l |Podesta and GaTy] l2tJTTl ISalem et al\ I2QT2). Ob- 
servations indicate that the turbulence is anisotropic and dominated by large 
perpendicular wavevectors, which is typical for KAWs possessing parallel electric 
fields and hence experiencing Cherenkov resonances. In the presence of a wide 
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KAW spectrum with overlapping harmonics, the protons undergo numerous 
Cherenkov kicks resulting in a Fokkcr-Planck proton diffusion in the velocity 
space along the mean magnetic field. 

We derived a quasilincar diffusion equation with the diffusion coefficient ex- 
pressed in terms of measured turbulence parameters: spectral break wavenumber 
/ii , turbulence amplitude W^i at the spectral break, and spectral slopes p2 below 
the break. At small and large /z, the asymptotic values for the spectral index are 
p ~ 1.7 and 2.7, respectively fe.o. lAlexandrova et aU I2TJ01?)) . There is however an 
intermediate spectral range, iii < /i < (12, where the spectrum slope is variable 
and often steepens to p ~ 3 — 4 or even higher values ( |Sahraoui et al., 2010[ ). The 
turbulence power at these scales is mostly in form of weakly/mildy dispersive 
KAWs that are Cherenkov-resonant with protons in the velocity range covered 
by observed nonthermal tails, V z cz (1 — 3) Va- Enhanced nonlinear interaction 
among KAWs can explain such steep spectra in this range (Voitenko and De 
Keyser, 2011), but resonant generation of nonthermal tails can contribute to 
the steepening as well. On the other hand, an additional source for the KAW 
replenishing in the dissipation range can be provided by the non-local nonlinear 
coupling among KAWs and large-scale MHD AWs QZhao, Wu, and Lu, 2011[ ). 

IRudakov et al.\ (|2012p applied an asymptotic analysis to study the quasilinear 
proton diffusion in the strongly dispersive KAW turbulence with /i> 1, which 
implies relatively high-energy resonant protons with V z 3> Va. However, such 
asymptotical analysis is inapplicable to typical tails with proton velocities Va < 
V z < 3Va- The reason is that the KAW phase velocity (Equation ©) cannot be 
approximated by a simple power-law dependence in the resonant wavenumber 
range 0.3 < \i < 3. We focused on this tail-resonant range, where most of 
kinetic-scale power is residing and where the full KAW dispersion (Equatiion ©) 
has to be used. Consequently, the resulting diffusion equation is more complex 
and does not have an immediate analytical solution similar to that found by 
IRudakov et al.\ (|2012{) . To proceed further analytically, we took into account a 
fast decrease of the diffusion coefficient with growing V z (see Figure 2), which 
allowed us to formulate a simple yet plausible model for the quasilinear plateau, 
Equations (|21[) and (pH]) , with the time-dependent front velocity V max and height 
/pi. By the use of this model in the diffusion equation (Equations (fT5|) and (flT)|0 . 
we found an analytical solution (Equation (|24|) ) for the plateau formation time 
as function of the plateau front velocity V max . Several numerical tests (one of 
them is shown in Figure 3) have confirmed the applicability of this model in the 
solar wind conditions. 

In general, our analysis and numerical estimations have revealed that the 
quasilinear proton diffusion driven by the observed solar wind turbulence pro- 
vides a robust generation mechanism for nonthermal proton tails, which can 
explain their routine observations in the solar wind. By implementing the mea- 
sured characteristics of the turbulence at proton kinetic scales fix, W^i, and pi 
in our model, we estimated the time scale of the quasilinear plateau formation 
in the velocity space. For typical nonthermal tails with V max = (1.5 — 2.5) Va 
the estimated formation time is 0.5 — 3 h, which is much less than the solar wind 
expansion time at 1 AU. This indicates that the local generation of such tails 
observed in the solar wind is a natural consequence of the observed kinetic-scale 
Alfvcn turbulence. 
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Some longer tails may require more time to develop, in which case a non-local 
radially-dependent problem has to be solved numerically using suitable solar 
wind models. In principle, with sufficient KAW power at /i > /Z2, the quasilincar 
plateau can extend well above 3Va, where an asymptotic /i » 1 analysis by 
IRudakov et al.\ (J2012I) is applicable. However, the turbulence level is usually 
insufficient to form the tails with V m!ix > 3Va within required time scales, which 
is reflected in the fast growth of the formation time beyond isw- This constrain 
is compatible with rare observations of such long tails. 

The most favorable conditions for the tail generation in terms of Vt p /Va occur 
in the range Vt p /Va — 1 — 1-5 for spectral slopes ranging from p2 =2.5 to 4, 
correspondingly. The optimal (Vt p /Va) t is not much varying with the spectral 
slopes in the range of interest, and we can say that in general (Vt p /Va) o . s=s 1. 
Such conditions are often met at 1 AU. As is seen from Figure 5, the range of 
Vt p /Va allowing for the tail generation has two bounds, upper (Vt p /Va) r2 and 
lower (Vt p /Va) ctV The upper bound (Vt p /Va) c1 . 2 is relatively large, its values for 
different spectral slopes are very scattered, and it is not very restrictive for the 
tail generation. On the contrary, the lower bound (Vt p /Va) ci1 is more restrictive. 
It has more collimated values that are not much different from each other and 
are not so small, which is seen from Figure 5. Say, for V max ~ 1-5Va tails 
(VT P /VA) crl — 0.25 and such tails cannot be generated locally if Vt p /Va < 0.25. 
The range of favorable Vt p /Va can not be extended towards smaller values. It 
seems that only short tails with V max < 1.5Va can be generated at Vt p /Va ~ 0.1. 
Relatively long tails are expected with steep spectral slopes P2 that are positively 
correlated with turbulence levels. 

We also analyzed a particular case reported by Sa hraoui et a l. (2010). We 
found that the measured turbulence can generate a nonthermal tail extending 
to Knax ;$ 2Va, but the proton number density in the tail is insufficient to 
explain the measured temperature anisotropy of the protons. We suggest that 
the nonlinear resonance broadening and enhanced collisional diffusion at V z < Va 
are important mechanisms significantly increasing the proton number density in 
the tail. Yet another possibility is that the longer tail producing larger measured 
anisotropy was generated at shorter radial distances. The related effects require 
further investigation. 

Behind the amplitude-related constrains mentioned above, the tail length 
and shape can also be constrained by plasma instabilities. Several instabil- 
ities can develop on the plateau-like proton tail, including the fast magne- 
tosonic, parallel Alfven-cyclotron (Rudakov et al., 2012[), and oblique Alfvcn 



( Daughton and Gary, 1999|| Voitenko and Goossens, 2003 ) instabilities. Although^ 
|Daughton and Gary|(|1999[) and IVoitenko and Goossensl (|2003p studied instabil- 
ity driven by the bump-on-tail, its generation mechanism is related to the anoma- 
lous Doppler (proton-cyclotron) resonance rather than the inverse Landau damp- 
ing, and hence is not restricted to the bump-on-tail but can also be efficient with 
plateau-like distributions. The cyclotron diffusion of tail protons, driven by these 
instabilities, can constrain tail lengths and produce proton-cyclotron quasilinear 
plateau reported bv lMarsch and Tul(|2001[) and lMarsch and Bourouainel(|201ip . 
or even produce the bump-on-tail features. In turn, the nonthermal proton tails 
and beams can re-emit KAWs ( |Daughton and Gary, 1 999} | Voitenko and Goossens, 2003| [ 
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|Wu, Chen, and Wu, 2012||Nariyuki, Hada, and Tsubouchi, 2012[ ) and significantly! 
modify the nonlinear evolution of the high-amplitude circularly-polarised Alfven 
waves (|Nariyuki, Hada, and Tsubouchi I2009I20I2"]) . A complex interplay of these| 



processes requires further investigations. 

Obtained analytical expressions for the proton diffusion coefficient can be 
incorporated in more sophisticated solar wind models accounting for external 
forces (electric and gravity), and both the particle-particle (Coulomb) and the 
wave-particle (Cherenkov) interaction terms in the right-hand side of Equation 
(|9"|). Because of their complexity, such models need to be solved numerically. 
Our preliminary simulation results demonstrate that the kinetic-scale Alfvcnic 
turbulence can create such nonthermal tails in the proton VDFs at the radial 
distances < 20 solar radii (see detailed descriptions of these simulations in our 
accompanying paper ( |Pierrard and Voitenko, 2012] )). 

We did not concern here the bump-on-tail features (proton beams) often 
observed in the solar wind. Such distributions have not been reproduced by 
the uniform quasilincar diffusion we studied. However, with variable level of 
the turbulence launched through the solar wind base, the irregularity of the 
generated tail length can easily produce the bump-on-tail features by the time- 
of-flight effect. Namely, in a region of enhanced turbulence the high-velocity 
end of a dense nonthermal tail propagates faster than the bulk protons. As 
a consequence, at larger radial distances this proton population penetrates a 
more quiet region with a weaker tail and appears there as a bump on tail. This 
mechanism is compatible with Helios observations of time-variable bumps on 
virtually persistent tails ( |Marsch et al, 1982[ ). 

In conclusion, appearance of nonthermal proton tails and beams in the solar 
wind is unavoidable in the presence of kinetic-scale Alfvcnic turbulence with 
sufficient amplitude. There are several other generation mechanisms for proton 
beams and tails (see Introduction), which can compete with each other or replace 
in varying solar wind conditions. More detailed statistical correlation studies are 
needed to discriminate among them and to find what mechanism is dominant in 
the solar wind. It would be particularly interesting to have a look at correlations 
between the tail lengths and/or densities and the turbulence levels at kinetic 
scales. 
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